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Epitaxial self-assembled quantum dots (SAQDs) are of both technological and fundamental interest, 
but their reliable manufacture still presents a technical challenge. To better understand the forma- 
tion, morphology and ordering of epitaxial self-assembled quantum dots (SAQDs), it is essential to 
have an accurate model that can aid further experiments and predict the trends in SAQD forma- 
tion. SAQDs form because of the destabilizing effect of elastic mismatch strain, but most analytic 
models and some numerical models of SAQD formation either assume an elastically homogeneous 
anisotropic film-substrate system or assume an elastically heterogeneous isotropic system. In this 
work, we perform the full film-substrate elastic calculation. Then we incorporate the elasticity cal- 
culation into a stochastic linear growth model. We find that using homogeneous elasticity can cause 
errors in the elastic energy density as large as 26%, and for typical modeling parameters lead to 
errors of about 11% in the estimated value of average dot spacing. We also quantify the effect of 
elastic heterogeneity on the order estimates of SAQDs and confirm previous finding on the possibility 
of order enhancement by growing a film near the critical film height. 
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I. INTRODUCTION 

Self-assembled quantum dots (SAQDs) func- 
tion as artificial atoms embedded in a semicon- 
ductor matrix.^ As such, they are useful for a 
rang e of electronic and optoelectr onic applica- 
tionG .'^'^"^''^'^"^''^'^'^'^"'"'^^'^^'^^'^^'^^'^^'^^'^^' For this 
reason, there has been a great deal of mod- 
elin g work on their dynamic formation pro- 

^,^g j20l21l22l23l24l25l26l27l28l29l30'l31l32l33l34l35l36l37l38l39l40l41l 

SAQDs are fabricated by depositing a semiconductor 
film on a lattice mismatched substrate with a smaller 
band gap, the most well-known examples being Ge^Sii-^: 
deposited on Si and lUj-Gai.^jAs deposited on GaAs. A 
good quantitative model of the SAQD formation can help 
in aiding understanding of the SAQD formation process 
and enable a sophisticated quantitative interpretation 
of experimental data, but more importantly, it can help 
move modeling from a descriptive mode to a predictive 
mode that could be used for process design optimization 
to aid in tasks such as the formation of new structures, 
control of morphology and enhancing order and repro- 
ducibility. Here, we imp rove upon previous spect ral 
models of SAQD formatioI^2a^^ ^^ ■ ^3 | 24 | 26 | 35 | 36 | 3 Il^ by 

incorporating elastic anisotropy and elastic heterogeneity 
simultaneously. We also estimate the errors introduced 
by neglecting these effects. 

Many reports in the literature make app roximations 
such as assuming elastic isotropyj ^" * ^ ^^'^^^^^ elastic h o- 
mogeneity of the film-substrate systenJ2iM1 36 | 37 | 39 | 40 | 



making a thin film approximation .^221231^32.3^ -^^g present a 
linear stochastic model of SAQD formation that incorpo- 
rates anisotropic elasticity and the elastic heterogeneity 
of the film-substrate system. We investigate the SAQD 
spacing (a mean property) as well as the order of SAQD 
arrays that is determined by the fluctuations in the spac- 
ing and alignments (Fig.[T]). Furthermore, we investigate 
the amount of error that previous approximations make 
(Table |l]). 

While there are other aspects of SAQD modeling that 
can be improved or incorporated, the presented work is 
an indispensable step in moving toward a more quan- 
titatively accurate SAQD formation model. The elas- 
ticity portion of the calculation applies generally to dif- 
ferent material systems, but other parts of the calcula- 
tion such as surface energies and diffusional dynamics are 
specific to group IV elements that have four-fold sym- 
metric SAQD formation dynamics such as Ge^^Sii-jj/Si, 
|36[37 43 g^j-^j j-^Q^ II- VI systems or III-V systems such as 
Ina;Gai_j;As/GaAs. 

SAQDs result from a transition from two-dimensional 
film growth to three-dimensional growth in strained epi- 
taxial films. When a flat strained film is perturbed by a 
film height undulation, elastic energy is released. When 
the released energy is greater than the cost in surface 
and wetting energy, the perturbation grows. This phe- 
nomenon is kn own as the ATG (Asaro-Tiller-Grinfeld) 
instability.'^^'^ Eventually the surface perturbations ma- 
ture into 3 D quantum dots. At a later stage the dots 
ripen, '2SESI although theoretically, they might form a uni- 
form array under some circumstances.^ ^^ * '^^ * ^^ * '^^ * ^^ * ^^ ! 
The interplay of the elastic energy, the surface energy 
and the wetting energy determines the energy landscape 
that drives SAQD formation, and the spectral modeling 
method used here yields a very transparent description 
of this interplay. 

The spectral model can be used to define and estimate 
parameters characterizing SAQD formation such as the 
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FIG. 1: Initial formation of Ge/Si SAQDs (Sec. |TT| with ex- 
aggerated height fluctuations for clarity. Lq is mean dot spac- 
ing. 2Lcor is the length over which SAQDs appear periodic. 
(Sec. 



nil 



characteristic length and time scales, the mean SAQD 
spacing, the alignment of SAQDs in an a rray as well as 
the critical film height for SAQD formation .I^SEIMIIIEIEI] 
In the stochastic form, spectral modeling can also elu- 
cidate order and reproducibility of SAQD arrays.!^ In 
higher order versions of spectral models known as a 
multiscale-multitime analyses, they can even elucidate 
longer term evolution of SAQDs. "^^ The clearly defined 
parameters from these models also inform finite element 
based model^^ and provide bench-marking for their per- 
formance. 

Elasticity is the most well understood influence on 
SAQD formation. As such, making fewer approxima- 
tions about elasticity will help investigations into other 
influences on SAQD formation that are more difficult to 
unde rstand. For example, wetting energy is barely under- 
g^QQ,j30i5015iJ^ and surface energy is generally treated as a 
constant, even though it almost certainly has strain and 
temperature dependence. It is also controversial as to 
whether surfaces should be treated as facets or not.'^^Ell 

While the importance of getting accurate estimates 
for a quantity as basic as the mean dot spacing is self- 
evident, the significance of SAQD order deserves further 
discussion. The ordering of SAQDs has been matter of 
concern in fostering the development of quantum dot 
based devices.^ There are two types of order, spatial and 
size. Spatial order is concerned with the uniformity of the 
spacings between the SAQDs and size order is concerned 
with the uniformity in the size of the SAQDs. The size 
and spacings of these self-assembled quantum dots are 
related, as dot volume is limited by the locally available 
material. Understanding what factors affect the order of 
SAQDs can guide experiments and simulation efforts and 
help in interpreting experimental and simulation results. 



Our enhanced elasticity calculation improves recent mod- 
els of SAQD order. We defer pattern fidelity in di- 
rected self-assembly to later work alt hough some initial 
results have been previously reported.'^SEIl 

Previous modeling work has focused on how elastic 
anis otropy and elastic h eterogeneity affect SAQD forma- 
^i^^mmmmmmmmb^t the two influences have been 
treated separately. The effect of elastic anisotropy has 
been studied in great detail. In Ref.f^^ it was shown that 
for heteroepitaxial system such as Sii-jjGe^^/Si, the sur- 
face undulations are likely to grow in the (100) direc- 
tions. It was also shown that for anisotropic materials 
the growth rate of the amplitude of the surface fluc- 
tuations is maximum when the wavelength is 4/3 the 
cutoff wavelength, similar to the isotropic approxima- 
tion. However, in the presence of a strong wetting effect, 
this ratio increases to 2. In the absence of misfit dis- 
locations, the islands are aligned in the (100) directions. 
However, experiments reveal that for films with thickness 
greater than the critical thickness for dislocation forma- 
tion, in the later stages of island formation, the islands 
align along the (110) directions due to formation of mis- 
fit dislocations.'^ In Ref.l^ a numerical investigation was 
carried out to study the effect of anisotropic strength 
on the formation, alignment and average isl and spa cing. 
More recent analytic studies on SAQD ordep^SEIISl com- 
plement these numerical studies. The effect of elastic het- 
erogeneity, however, has received more limited attention. 
In Ref.42j a linear stability analysis was performed that 
incorporated the elastic stiffness for both film and the 
substrate. One major conclusion was that elastically stiff 
substrate has stabilizing effects on the film that dimin- 
ishes with increasing film thickness. In Ref.^^ a nonlinear 
evolution equation was derived using a thin-film approx- 
imation. However, Refs.l^ ancP^ approximate elasticity 
as isotropic. 

Here, we treat elastic effects without approximations 
regarding isotropy, homogeneity or film thickness. We 
find various parameters that can be derived and esti- 
mated from spectral SAQD growth models, and we com- 
pare them with the results of the other more approximate 
models (Table |l]). These parameters are £, the elastic 
energy density coefficient (Figs. [2] and [3| Le, the per- 
turbation wavelength that is the most energetically un- 
stable (Fig. [4]), Lqj the perturbation wavelength that is 
kinetically most unstable and gives the mean dot spacing 
(Fig. [6|, and Ucor (Fig. [?]), the number of dots in a row 
whose positions are well correlated. Each of these values 
is compared with the predictions of more approximate 
models, namely, the elastically anisotropic homogeneous 
approximation, the elastically anisotropic thin-film ap- 
proximation, the elastically isotropic heterogeneous ap- 
proximation, the elastically isotropic homogeneous ap- 
proximation and the elastically isotropic thin-film ap- 
proximation presented recently.!^ For the order analy- 
sis, ncor, comparisons are only made with the elastically 
anisotropic models as elastically isotropic models are not 
suitable for order predictions of periodic arrays.^ Also, 
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all of the reported estimates depend on the average film 
height (H); thus for each comparison we present a calcu- 
lation corresponding to a typical average film height of 
H — 4.25 ML — 1.2 nm in Table m with some additional 
values displayed in Figs.|3] [4] and [6]-[7] It is worth noting 
that all of these approximations correspond to various 
limits of our elasticity calculation. For example, the ho- 
mogeneous anisotropic approximation is identical to the 
limit as the average film height {H) becomes large. The 
anisotropic thin-film approximation corresponds to the 
limit as _ff — !■ 0, and the various isotropic approximations 
can be obtained by using an isotropic elastic stiffness ten- 
sor by, for example, taking the Voigt or Reuss average of 
the actual elastic moduli. Finally, we give in-depth anal- 
ysis throughout only for the anisotropic models. 

We model the SAQD growth process with a stochas- 
tic surface diffusion model. We perform a linear analysis 
of the dynamics of the film evolution, which corresponds 
to small height fluctuations. Although such an analysis 
would only be valid for the onset of island formation, it 
determines the initial placement of SAQDs; thus deter- 
mining the initial mean-spacing (Lq) and ord er (ricor)- At 
later stages, the SAQDs either order or ripeipIlSMIIllIlll. 
The spacing and order established at the small fluctua- 
tion stage will influence the order at a later stage. This 
has also been verified through non-linear calculations in 
Ref.'SSl. Linear effects also set the length and the time 
scale for measuring the pcrturbations,SZl and determine 
the arrangement of dots. Linearization offers a transpar- 
ent way for analysis and is also a prerequisite for under- 
standing more advanced non-linea r model s. The proce- 
dure for order analysis follows Refs.'^^'^^'^. Most models 
in literature are deterministic; however, the stochastic 
model is more realistic, as there is no rigorous physical 
explanation for the artificial initial random roughness in 
the deterministic models. 

The rest of this article is organized as follows. We 
give details of the stochastic surface diffusion model in 
section ini In section ITlTl we discuss the order calculations 
using film height correlation functions. We present our 
conclusions in section HVl 

II. MODEL 

The formation of SAQDs takes place through surface 
diffusion that is driven by a diffusion potential /i and 
contains thermal fluctuations, ^(x, t)!^. /j, is a non-local 
functional of the film height H and a function of the hor- 
izontal position X = {x, y) (Fig.[T]), so that /i — > ^jl[H] (x). 

The normal velocity of the evolving film surface is 

v^ = DVl^i + Vs■iM, (1) 

where is the surface laplacian, Vs- is the surface di- 
vergence, and we omit explicit coordinate and time de- 
pendences for brevity. Here we consider the case of an- 
nealing of a film and therefore we omit a surface flux term 
in Eq. [l] 



We linearize all quantities about the average film 
height H , 

H{^,t) = H + h{ii,t), (2) 

where the average film height H can be controlled by con- 
trolling the amount of deposited material and /i(x, t) rep- 
resents the fluctuations about this average that cannot be 
experimentally controlled. In this procedure, the elastic 
contribution is non-local, so analysis is ai ded b y working 
with Fourier components. Following Refs.l^^^^we use the 
Fourier transform convention, /(x) = / (Pxe^^'^fk and 
/k = (27r)~^/ (fke-'^''f{x). Note that subscript k is 
used to indicate functions of wave vector k while (x) is 
used to indicate dependence on the real-space coordinate. 
We proceed in two steps. First, we linearize the diffusion 
potential ii. Then, we linearize the dynamic governing 
equations. Similar to Ref.'^, we keep terms only to lin- 
ear order in /i(x,i). 

Previously, the homogenous elasticity approximation 
was used t o ide ntify three related wavenumbers and 
wavelengthsj ^'* * ^^ ^ the characteristic or cutoff wavenumber 
and wavelength kc and Lc = 27r/fcc, the wavenumber and 
wavelength for maximu m energy release, — (1/2) fc^ 
and Le — 2Lc^^^^^ and the wavenumber and wave- 
length of th e fast est growing mode was identified, fco ^^nd 
Lq — 27r/fco.'2311! In the absence of a wetting effect such 
as for thick films, ko = (3/4) kc {Lq = i/SLc), while for 
thin films where the wetting effect is strong, kg ranges 
from ko — kE {Lq — Le) at the crit ical film height to 
ko = (4/3)/c£; {Lo = (3/4)L£;).l3Il33l3lin the less approx- 
imate formulation that is elastically heterogeneous and 
anisotropic, these relationships are not as simple. In the 
following analysis, we identify fc^ and ko- 
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FIG. 2: The elastic energy density prefactor as a function of 
the dimensionless vector for Ge/Si at 600 K. 
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A. Energetics 



The difTusion potenti al, fi, consis ts of three parts, 
M = /iciast. + Msurf + Mwet. ■ '^^'^^ ' ^'^ ' '^ The elastlc energy 
part destabihzes the two-dimensional growth made, the 
surface energy term stabilizes the short wavelength (high- 
k) modes, and the wetting potential stabilizes all wave- 
lengths. We proceed by calculating the Fourier transform 
of the diffusion potential, /ik, to linear order in terms of 
the Fourier transform of the film height, /ik- 



1. Elastic amsotropy and heterogeneity 

The elastic contribution to the diffusion potential is 
just the elastic energy density at the film surface, de- 
noted uj (x) times the atomic volume (/ieiast. = i^Lu{^))^^ 
We proceed by calculating the Fourier transform of w(x), 
ciJk to linear order in surface height fluctuations, /ik while 
taking into account the effect of both elastic heterogene- 
ity and elastic anisotropy. 

The full calculation is described in the Appendix, and 
it results in an elastic energy density of the form 



(3) 



where £g^^{kH) is the elastic energy density prefactor that 
depends on both wave vector direction ^k E^nd dimen- 
sionless product kH . Note that this form is equivalent 
to writing £ (ki?) where kTJ is a dimensionless vector 
quantity. This result should be contrasted with previous 
calculations. In the homogeneous isotropic approxima- 
tion, the prefactor is a constant, and in the homogeneous 
anisotropic approximation, th e pre factor depends only 
on the wavevector direction g^. p4 | 37| 

We perform numerical calculations for [001]-oriented 
Ge/Si at 600 K to give a concrete example of the en- 
ergy prefactor f ('c^) • The elastic stiffness tensor djki 
is 4 — fold symmetric for rotations about the [001] axis; 
thus £g^^ (^^) is also 4 — fold symmetric. This symme- 
try manifests itself in the arrangement o f SAQDs into 
a four-fold symmetric quasiperiodic latticd^lIllIlTMl^ ■\Ye 
use the following physical constants. The superscripts / 
and s differentiate between the elastic constants of the 
film and the substrate respectively. For Ge at 600 K, the 



elastic constants are c{i = 11.99 x lO^^dyn/cm^ 
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4.01 X 10"dyn/cm2 and = 6.73 x 10" dyn/cm^.l^s! 
For Si at 600 K, cj^ = 15.61 x lO^dyn/cm^, cl^ = 
5.63 X 10"dyn/cm2 and 0%^ = 7.82 x 10" dyn/cm^.ISi 
Using cgc — 0.5658 nm and agi = 0.5431 nm, the mis- 
match strain is = 0.0418. Figure [2] shows a plot of 
the elastic energy prefactor, £g^^ {^^^ against the dimen- 
sionless variables fca;ff and fc^ 5. Figurejsjshows fg^ {^^) 
as a function of kH for three values of 0k along with a 
comparison to the discussed homogenous and isotropic 
approximations. In Figure [3] we can clearly see that 
as kH — > oo, the prefactor reaches its asymptotic values 
that also correspond to the homogeneous approximation. 
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FIG. 3: Elastic energy density prefactor fg^ for Ge/Si at 
600 K from isotropic approximation^^ and anisotropic cal- 
culations for &k = 0°, 22.5° and 45°. Asymptotically large 
kH limits corresponds to the anisotropic/isotropic homoge- 
neous approximations and are shown as dotted lines (...). 
Percent error in the values of £"9^ for = 0° is shown for the 
anisotropic homogeneous approximation. 



Typical values for film height are H < 20 ML, while 
typical relevant wavelengths are 30 — 40 nm; thus, rele- 
vant values for kH are < 1 . In figure [3] we also plot the 
error due to the homogeneous anisotropic approximation 
for ^k = 0° and < kH < 2.5. For values of kH > 2.5 
the error is significantly less (error < 1%). It should be 
noted that we focus primarily on the values at 6'k = 0° 
because undulations are more likely to grow in the (100) 
directions^^^^ At lower kH values {kH < 0.4) the er- 
ror is higher (> 10%) with the upper bound being 26%. 
For example, for a periodic array of islands spaced at 
Lq = 49.06 nm(fc — 2tt/Lq — 0.128 nm~^) with average 
film height H = 4.25 ML = 1.2 nm so that kH = 0.154, 
the error in the calculation of elastic energy density is 
about 18%. We find that the anisotropic thin-film ap- 
proximation does a bit better with an error of —7%. 
We report these final values along with comparisons to 
other approximations in Table. [T] Such errors limit the 
accuracy of quantitative models, as this error propagates 
to calculations of the various characteristic wavenumbers 



and critical film 



(Sees. II A 3 and II B I, mean dot spacing, rate of growth 



height. 



2. Surface and Wetting energies 

The other contributions to the SAQD formation ener- 
getics are the surface and wetting energies. Since our fo- 
cus is on 4-fold symmetric systems, the only anisotropic 
term is due to the elastic energy.l^II As in Ref.l23, 



Msurf.,k = ri(7fc^)ft.k, 



(4) 
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TABLE I: Comparison of presented model with various approximations. Our mo del us es heterogeneous amsotropic elasticity 
(Het. Anis.). Other models use homogeneous anisotropic elasticity (Hom. Anis.j^^'^, anisotropic thin-film elasticity (Anis. 
Thin), heterogeneous isotropic elasticity (Het. Iso.}^- , homogeneous isotropic elasticity (Hom. Iso.) and isotropic thm-film 
elasticity (Iso. Thinj^ All calculations use average film height H = 4.25 ML and fo° (0.154) uses k = 27r/(49.06nm) = 
0.128 nm~^. Values in parentheses indicate % error due to each approximation. 
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(% error) 


(% error) 


(% error) 


(% error) 


(% error) 


Eo" (0.154) 


4.18 


4.92 


3.90 


5.64 


6.47 


5.24 


(10^ erg/cm^) 




(+18%) 


(-7%) 


(+35%) 


(+56%) 


(+25%) 


Le (nm) 


55.4 


49.2 


62.1 


39.9 


37.4 


46.2 






(-11%) 


(+12%) 


(-28%) 


(-32.5%) 


(-17%) 


Lo (nm) 


49.1 


43.7 


55.1 


35.4 


33.2 


41.0 






(-11%) 


(+12%) 


(-28%) 


(-32%) 


(-17%) 


'^cor. 


2.71 


2.07 

(+24%) 


2.62 

(-3%) 


n/a 


n/a 


n/a 



whe re 7 can be interpreted as the effective surface en- 
ergy.l2Z13fl ^pj-^e linearized wetting potential is 



(5) 



where W" is the second derivative of the wetting poten- 
tial with respect to the film height evaluated at the aver- 
age film height H = H. For the example here, we follow 
Ref.'2il and take the wetting potential to he W = B/H, 
where B is a material constant. 



3. Energy cost function 

Combining Eqs.[3]|4]and[5] we can write the linearized 
diffusion potential in Fourier space as 



= /(fc, 6'k,H)/lk, 



(6) 



where /(fc, 9^, H) ^ n [-k£g^{kH) + -fk^ + W"] is the 
energy cost per unit height for a periodic perturbation. 
The minima in the energy cost function lie along the 
(100) directions (^k = 0°, 90°, 180°, 270°) and occur 
at wavenumber ks so that the most energy is released 
when dots form at a period of Le — 2TT/kE- Le is 
a function of H, a dependence that is due purely to 
the more precise elasticity calculation we present, and 
not, for example, a result of the wetting potential. For 
a concrete example, we use the estimated surface en- 
ergy density 7 = 1927 erg/cm^ and the atomic volume 
= 2.27 X 10"^'^cm'^. Figure |4] shows Le as a function 
of H along with its values for the anisotropic homoge- 
nous and isotropic approximations. As H becomes large, 
Le approaches the anisotropic homogeneous approxima- 
tion value that is independent of H. For an average film 
height, H = 4.25 ML = 1.2 nm, the error in Le from the 
homogeneous approximation is about 11%. We report 
results for all approximations in Table. [l] 




FIG. 4: Characteristic wavelength Le as function of H for 
Ge/Si at 600 K from isotropic and amsotropic calculations. 
Corresponding homogeneous approximations are shown as 
dashed lines. The percent error in is shown for the amsotropic 
homogeneous approximation. 



The energy cost function is also useful for determin- 
ing the critical film thickness for SAQD formation, or for 
modeling purposes it can be used to estimate the wetting 
potential that leads to an observed critical height. The 
critical height for SAQD formation in the Ge/Si film- 
substrate system is generally observed to be 4-6 ML.ISSl 
Here, we choose a critical film height of 4 ML and follow 
the procedure from Ref.'^. We assume a wetting poten- 
tial of the form W{H) — B/H and then find the coeffi- 
cient B that gives a critical film height. He — 4 ML, by 
setting the minimum value for the energy cost function to 
zero, /mill — /{ks, 0°, i/c) — 0. Solving for B is a simple 
procedure as / is linear in B. We find that even the value 
of B is sensitive to the anisotropic homogeneous and 
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other approximations. At 4 ML, B = 1.61 x 10~^ erg/cm 
for the full theory, and B = 2.28 x 10~^ erg/cm for the 
anisotropic homogeneous approximation, about 42% dif- 
ference. 



B. Dispersion Relation 

The linearized evolution equation in Fourier space is 
given by2Sl 



5t/ik = CTk/ik + V^^DhT [{k ■ r7k(i)] , (7) 

where the second term is the Fourier transform of ^(x, t) 
to linear order, {Vw{t)vi'it')) = {2n)-'^iS^{k-k')5{t-t') 



(8) 



is the generalized dispersion relation that gives the rate 
of growth (positive values) or decay (negative values) of 
each height Fourier component /ik. 



ak/D (lO^erg/cm 




ky (nm ^) 



FIG. 5: o-k/D vs k for // = 4.25 ML for Ge/Si at 600 K. 

Figure [s] shows the dispersion relation Ck for H = 

4.25 ML. For the case shown in figure [s] CTk has 4 peaks 

along the (100) directions at ko = (0, ±0.128) nm"! 

and (±0.128, 0) nm~^. The four peaks indicate that 

the instability is maximum in the (100) directions thus 

making them the likely directions for the alignment of 

S AQDs. Thi s alignment is consistent with previous stud- 
jggPmiinilTl 

Similar to Ref.l^^, we expand dk about its peak values 
to get 



O'O - ^cr||(^ - fc||)^ - ^C^-L^i, 



(9) 



where 



aVk 



dk^ 



1 aVk 



[ko,{eo),^] 



k^ 



06^ 



[koA0o) 



n corresponds to the number of peaks, Oq is the orienta- 
tion of kg, and k\\ and fc_L are components parallel and 
perpendicular to kg. We discuss the dependence of fep 
and the mean dot spacing on film thickness next along 
with the discussion of SAQD array order. 



III. ORDER ANALYSIS 

The spatial order of SAQDs is best characterized by 
the mean geometric spacings, Lq, and alignments and by 
the degree of fluctuation about these means. The aver- 
age alignment of SAQDs is (100), and we characterize 
the range of order by, neon the number of dots in a row 
whose positions are likely to be well correlated, mean- 
ing that they are likely to be both regularly spaced and 
well- aligned. In the following discussion we present cal- 
culations for different average film heights, and for each 
film height we calculate average dot spacing and the num- 
ber of correlated dots when film height fluctuations reach 
atomic scale size. For the second part, for finding ricor, 
we use the film height correlation function and ass ociate d 
correlation lengths which were derived previously.^^^ 



A. Average Dot Spacing 



As done previously,'2il2ll2lIlZIl2l estimate the aver- 
age initial spacing between dots to be Lq = 27r/fco. Fig- 
ure [g] shows a plot of Lq against H and compares it with 
the results for the anisotropic homogeneous and isotropic 
approximations. The error associated with the homo- 
geneous approximation is also shown. We report val- 
ues for H = 1.2 nm in Table |l] For the example stud- 
ied here (Ge/Si at 600 K), the value of average spacing 
for anisotropic heterogeneous elasticity calculation varies 
between 32.8 nm to 55.7 nm. Note that as -ff — > cx), 
Lq reaches its asymptotic value, which is the same as 
the value from the anisotropic homogeneous approxima- 
tion. Thus, the anisotropic homogeneous approximation 
is valid for thick films. Typically experiments correspond 
to values of H that are less than 20 ML(4.25 nm). 



B. Order Analysis Using Correlation Functions 

The autocorrelation function and its Fourier trans- 
form, also known as the po wer spect rum, are very useful 
for characterizing dot order .IS^EiMl xhe autocorrelation 
function is defined as 



C^(Ax) 



(i^x/i(x + Ax)/i(x). 



(10) 



For an imperfectly periodic array of SAQDs the auto- 
correlation function decays away from the origin. The 
distance over which the autocorrelation function decays 
is known as the correlation length Lcor- The value 2Lcor 
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limiting. Taking the inverse Fourier transform, the cor- 
relation function is 



C(Ax) 



DnhTP 



,2(70* 



■^{Ax'/Ll+Ay'/Ll) 



cos(A:Aa;) 



7r(Toi||ij_ 

+e-^(A-Vii+A«ViS)cos(/fcAy)l . (12) 



//(nm) 



Figure [7] shows the number of correlated dots calcu- 
lated as ricor = 2L_l/Lo for the small fluctuation stage 
(C(Ax = 0) = 1 ML^). Both error and number of cor- 
related dots decline sharply for a small increment in film 
height above the critical film height. We find the error 
drops from 24% at = 4.25 ML to_3% at = 4.95 ML. 
Note that with further increase in H, the error fluctuates 
between ±2% before reaching 0% for higher H values. 



FIG. 6: Average dot spacing Lq as function of H for Ge/Si 
at 600 K from isotropic and anisotropic calculations. Corre- 
sponding homogeneous approximations are shown as dashed 
lines. The percent error in is shown for the anisotropic homo- 
geneous approximation. 



represents the distance over which the SAQDs appear to 
be periodic meaning regularly spaced and well aligned. 
The power spectrum is 

The power spectrum for a nearly periodic array of SAQDs 
will have peaks with finite width, A^. The spectrum 
peak width A^ is related to the correlation length Lcor 
by Lcor = 1/Afe. 

Each simulation or experiment corresponds to one par- 
ticular realization with its own autocorrelation function; 
however, for sufficiently large simulation sizes, the fluctu- 
ations in (Ax) are small, and the ensemble average of 
the autocorrelation functions can be predicted and pro- 
vides a good estimate of individual autocorrelation func- 
tions and spectrum functions.'^ The ensemble average of 
the autocorrelation function is the correlation function 
C(Ax) = (C'^(Ax)) = {h{Ax)h{0)). Similarly, the en- 
semble average spectrum function is Ck = (C*^), where 
Ck is also the Fourier transform of C(Ax), and Ck is 
the coefficient in the covariance of the Fourier compo- 
nents ft-k; (^k^k') = Ck<5^(k — k'), where (5^(k — k') is the 
two-dimensional Dirac Delta function. 

T he spe ctrum function can be solved using Eqs. [7] [9] 
and [Tops 



DnhT _ 



n=l 



-i 



(11) 



where L|| = ^J2a^ and = \/2a±t are the correla- 
tion lengths. L|| gives about half the length over which 
the dot spacing is regular, while L±^ gives about half the 
length over which a row of dots is straight. Of the two 
correlation lengths, L±^ tends to be smaller and thus more 
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FIG. 7: Number of correlated dots ricor for small height fluc- 
tuations Ge/Si at BOOK for anisotropic heterogeneous and 
anisotropic homogeneous calculations. Also shown, error in 
ncor from anisotropic homogeneous approximation. 



IV. CONCLUSION 

Most theoretical and many numerical models of SAQD 
growth approximate film-substrate systems as elastically 
homogeneous. We have examined the effect of elastic 
heterogeneity on SAQD mean spacing and the order es- 
timates developed in Refs.l^ ancP^. We have performed 
a linear analysis incorporating both elastic heterogeneity 
and elastic anisotropy. We quantify the effect of hetero- 
geneity as percent error in the calculated values of elastic 
energy density, average spacing between the SAQDs and 
number of correlated dots based on homogeneous approx- 
imation. We show that the homogeneous approximation 
of the film-substrate system can lead to significant errors 
in the calculations for formation and ordering of SAQDs. 
For the case of Ge/Si system at 600 K the upper bound 
for error in the calculated value of elastic energy density 
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is found to be as large as 26%. The error declines as H in- 
creases. For a typical average film height, H = 4.25 ML, 
we calculate an error of about 11% in the estimation of 
average spacing between the SAQDs. Using a stochastic 
model, and the film height correlation functions, we find 
that the error in the estimated number of correlated dots 
declines quickly as H increases. The error in the esti- 
mated number of correlated dots drops from about 24% 
at ^ = 4.25 ML to about 3% &t H ^ 4.95 ML. For thin- 
ner films, the thin-film approximations can reduce this er- 
ror, but error still remains as the thin-film approximation 
actually overestimates the effect of elastic heterogeneity. 
In general we find that the most error is due to using an 
isotropic approximation. For the isotropic heterogeneous 
approximatiorpS the error in mean dot spacing remains 
more or less constant at 34% for values oi H < 20 ML. 
We did not report order predictions from isotropic mod- 
els, as they are inappropriate for order and alignment 
estimates."^ 

The interplay between elastic strain, surface energy 
and surface diffusion can be quite complicated. Errors 
introduced by the elasticity portion of models can con- 
found our ability to asses how well we model surface and 
wetting energy. Given the challenge of developing accu- 
rate models of surface and wetting energies, it is essen- 
tial that the elasticity part of the calculation be correct. 
Here we have used a linear model that includes elastic 
anisotropy and heterogeneity. This is an important step 
in the development of a complete non-linear stochastic 
model required for a more comprehensive quantitative 
analysis, for example incorporating surface energy and 
diffusive anisotropy. We have demonstrated that one 
must incorporate both elastic anisotropy and elastic het- 
erogeneity to avoid the introduction of significant errors 
in the calculation of the elastic energy density and mean 
dot spacing. This effect will have important consequences 
for single layer arrays grown on a flat [001] surface, but 
should also suggest the magnitude of errors introduced in 
modeling other SAQD morphologies such as multilayers 
and growth on patterned substrate. 
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APPENDIX: ELASTIC ENERGY COEFFICIENT 

The increase in elastic energy due to the addition of a 
small material volume at the surface is just the elastic en- 
ergy density at the surface.!^ We calculate the elastic en- 
ergy density using perturbation theory following Ref.'^, 
but here we take elastic heterogeneity into account. This 
calculation was performed previously but only with the 
approximation of isotropic elasticity.^^ Here we incor- 



porate both elastic heterogeneity and elastic anisotropy 
to calculate the elastic energy density at the film surface. 

We consider a flat fllm on a substrate. The lattice 
mismatch between the film and the substrate introduces 
a misfit strain e„i in the film and leads to a uniform stress 
distribution in the film given 





(T„i 





where cr„ 



Me™, and M is the biaxial modulus, 



M = 



c{i + 0(2 2(c{2) /c(i 



(13) 



(14) 



where c\i and 0(2 are elastic constants for the film. We 
perturb the film surface so the film height fluctuates as 



ikx 



(15) 



where the cartesian coordinate system is set up so that 
2; = or the x-y plane lies at the interface of the fllm- 
substrate system, and the a;— direction is aligned along 
k. Then, we calculate the elastic energy to flrst order in 
the perturbation amplitude, /iq. This calculation requires 
four steps. First we must flnd the surface normal vector 
n to first order in Hq. Then, we must find the admissi- 
ble equilibrium eigenmodes for the elastic displacement 
that have the same periodicity as the height perturba- 
tion. Then we must find the eigenmode coefficients from 
the surface boundary conditions and internal matching 
conditions (compatibility and equilibrium). Finally, we 
find the elastic energy density at the free surface to first 
order in /iq. The first step is simple, to the first order in 
the normal to the surface of the film is given by 



n(x) = —ikhoe '-^Bx + e 



(16) 



The remaining three steps follow. 

We find the elastic deformation eigenmodes that have 
the same periodicity as the perturbation, and that sat- 
isfy internal equilibrium; working with the displacement 
field u automatically satisfies compatibility away from 
internal interfaces. We first construct the rank 4 elastic 
stiffness tensor for both film and the substrate for an ar- 
bitrary passive rotation 0k in the x-y plane so that k can 
lie along any direction in the plane. 

3 

■i j'jA:,/ — 1 

(17) 

and 

3 

1 

(18) 

where we are using the superscript / and s for the film 
and the substrate respectively and i?(0k) is the passive 
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rotation matrix by an angle 9^ about the 2:— axis. To 
match boundary conditions in a later step, displacements 
must have the form 



(19) 
(20) 



where k and C arc unknown eigenvalues. The stress ten- 
sors in the film and the substrate are 



s,t=l 

3 



^* dxt ^ 



s.t=l ^ 



(21) 
(22) 



For the film, the elastic equilibrium equations are 



^ a^cU(^k)^"f = 0;n = 1...3 (23) 

q,s,t — l ^ 



(24) 



\t=i 



where 



q,s—l 

To obtain non-trivial solutions, we set the determinant 
of C/t(6lk,K) to zero. We thus obtain six eigenvalues of 
K denoted by with p = 1 ... 6. Each value of k = k^" 



is substituted back into Crt(6'k,K), and Eq. 24 is solved 
to find the corresponding eigenvectors Uj' . The displace- 
ment components for the film in terms of the unknown 
coefficients Ap are thus 



(26) 



where are the eigenvectors, and Bp are the unknown 
coefficients. 

We now find the nine unknown coefficients [Ap and 
Bp) using the traction-free boundary condition at the 
surface, and the internal matching conditions at the film- 
substrate interface, namely equilibrium and compatibil- 
ity. The traction on the surface of the film is 



(28) 



q=l 



where z = H + h{x). Substituting Eq. [2lJinto Eq. ( p8| , 

we get 



E 

q,s,t=l 



dxt 



(29) 



3 6 

= ieraho ^ ^ c{j^,tkApUf{iSti . . . 

1 p— 1 

... + KP(5t3)e'=("+'"''^)n, + {a„i)rqnq. 

Again, we substitute z — H + /i(x), and we keep only 
terms up to first order in Hq to get 



Tr 



3 6 



E E H.l(^k) + A«^C3^..3(^k)) A, 



s—1 p—1 



(30) 



xC/f, 



kK^H 



M5rl 



Since the traction on the film surface must be zero, we 
have 



3 6 



EE0c3-l(^k) + '^^cL3(ek))A, 



jjP^n^kH ^ ^^^^^ 



s—1 p—1 



(31) 

giving three equations for r — 1 ... 3. 

The force balance at the internal film-substrate inter- 
face requires 



where we assume that the perturbing elastic field dis- 
placement components are proportional to em and ho, 
and we put the pref actor i = in for convenience. 

We use the same procedure to find the eigenvalues and 
eigenvectors for the substrate displacements u^, where 
C^f{9k,() has the same form as Eq. 25 but using the 
substrate elastic constants, c^^^j. Six eigenvalues, are 
obtained; however, we assume that the substrate is a 



semi-infinite solid so the displacement field uf 



at 



z — ~oo. Thus, we only retain the three eigenvalues 
with Re [(P] > that satisfy this condition and discard 
the other three. We find the displacement components 
of the substrate, 



ihoemY.BpV/'e''('^+i''^\ (27) 
p=i 



0-3r = O'Sr 



(32) 



In terms of the unknown coefficients Ap and Bp, we can 
write Eq. ( 32 1 as 



3 6 



3 3 



E E (*'^3rsl + '^^^3^53) ^P^s — E E ('"^Srsl • ■ • 
s—1 p—1 s — 1 p—1 

■ ■.+Ccl^,)BpVP 

(33) 

for r ~ 1 ... 3. For the compatibility between the film 
and the substrate at the interface, the displacements of 
the film and the substrate must be equal, so that 



(34) 
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In terms of the unknowns, the compatibiHty equation can 
be written as 



p=l q=l 



(35) 



giving three equations i — 1 ... 3. We then calculate 
the nine coefficients, Ap{0\^,kH) with p — 1...6 and 
Bq{9i^,kH) with q — 1...3 using using Eqs. (31 1, (33) 
and (|35l 



Following Ref.E3, the we find the elastic energy at the 
film surface to first order in ho to be 



(36) 



where Uq is the energy of the unperturbed flat film, a 
constant. Using Eq. [26| 



U = Ua- £e^{kH)khoe"''' , where 



(37) 



£e^{kH) = Me,2„ Vv4p(0k,fci?)C/f(^k), (38) 



where we note that Ap will depend on (?k, k and H , and 
Uf will depend on 0k- By the principle of superposi- 
tion, we can use the elastic energy coefficient £0^ {^H) 
for sums of periodic perturbations as well. 
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